Analytical approach to chiral symmetry breaking in Minkowsky space 
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The mass gap equation for spontaneous chiral symmetry breaking is studied directly in Minkowsky 
space. In hadronic physics, spontaneous chiral symmetry breaking is crucial to generate a constituent 
mass for the quarks, and to produce the Partially Conserved Axial Current theorems, including a 
small mass for the pion. Here a class of finite kernels is used, expanded in Yukawa interactions. The 
Schwinger-Dyson equation is solved with an analytical approach. This improves the state of the art 
of solving the mass gap equation, which is usually solved with the equal-time approximation or with 
the Euclidean approximation. The mapping from the Euclidean space to the Minkowsky space is 
also illustrated. 



I. INTRODUCTION 

Here I solve directly in Minkowsky space the mass 
gap equation for Spontaneous Chiral Symmetry Break- 
ing (S^SB). S^SB was introduced by the original work 
of Nambu and Jona-Lasinio 0, and it is presently ac- 
cepted to occur in hadronic physics, where it generates a 
constituent mass for the quarks. S^SB also implies the 
Partially Conserved Axial Current theorems |2|, includ- 
ing a small mass for the pion. In the literature the mass 
gap equation is usually solved either in equal time Q 
or in Euclidean space [J] , in order to avoid the poles and 
complex quantities which are expected in full Minkowsky 
calculations. Recently the scientific community has been 
exploring different approaches to the Minkowsky space 
5]. An exact solution in Minkowsky space will test the 
quality of the aproximate solutions. Moreover a solution 
in Minkowsky space opens wider applications. For in- 
stance the Bethe Salpcter amplitudes can be computed 
on the mass shell momentum p 2 = M 2 > 0, and at the 
same token it is possible to boost the hadrons to any 
convenient frame. In this paper an analytical approach 
is applied to finite and analytic kernel. For simplicity the 
infrared or ultraviolet divergences are not addressed here. 
Although a complete solution can only be achieved nu- 
merically, an interesting insight is nevertheless presented 
by the analytical approach. 

In the Schwinger-Dyson formalism truncated at the 
rainbow level, the mass gap equation is, 



S~ l (p) = S Q \p) - 
S(p) 



-iV(p-q) ■■ 
f S(q) ^ 7 



A(p 2 ) jp - B(p 2 ) + ie 



(1) 



For the quark-quark interaction I consider the class of 
kernels that can be decomposed in a class of Yukawa 
potentials, which is finite both in the ultraviolet and in- 
frared limits. 



-iV(p-q) = -iJ2 ^_ „S^2 



(p - q) 2 - X 2 + ie 



(2) 



To ensure that the kernels are finite in the ultraviolet 
limit of large (p — q) 2 I assume that on — 0, since this 
implies that the kernel vanishes at least like \(p — g) 2 ] -2 . 
In the infrared limit the kernels are finite when the Xi are 
finite, or when infrared cancellations occur. The class of 
kernels defined in eq. @ is quite general, for instance it 
includes the Coulomb interaction regularized by a Pauli- 
Villars term. It also includes the Fourier transform of 
the confining linear potential, 



F.T. 



(|r|e 



-in 



d\ 2 k 2 + A 2 



(3) 



and of the linear potential with a negative infinite shift 

d 4-7T 



F.T. 



dX k 2 + A 2 



(4) 



which are both obtained in the limit of vanishing A. Using 
finite differences, the derivatives in A can be decomposed 
in the class of potentials of eq. J5J. The class of kernels 
addressed in this paper is not only general, it is also co- 
variant, analytical and forward propagating, and this is 
convenient to study the mass gap equation in Minkowsky 
space. 

Chiral symmetry breaking occurs when a non- 
vanishing mass is dynamically generated. For clarity let 
us also assume that A ~ 1 in eq. Q). This approxima- 
tion is qualitatively acceptable when the kernel is finite 
0,0. Computing A is not difficult, but it obscures the 
result of the paper. Then the mass gap equation is a 
single non- linear and integral equation for B = M, 



M{p 2 



mo 



i d q 

-co Wj* 



V(p - q) 



M(q 2 



q 2 - M 2 (q 2 ) + it 



where all Dirac and colour algebraic factors are absorbed 
in the coupling constants ai of the potential V(p— q). The 
chiral limit of vanishing current quark mass (the mass in 
the free quark propagator) mo — 0, is also assumed. This 
is particularly interesting because it applies to the physics 
of the quarks u and d. It is clear that eq. (JSJ then has a 
trivial solution M(p 2 ) = 0. The problem that this paper 
addresses is the other possible solutions of eq. ©, with 
the kernel of eq. © • The technical difficulties reside in 
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the multiple integral with poles and complex quantities 
and in the non-linearity of the self-consistent equation. 

This paper is organised in Sections. In Section [H] 
standard approximations are applied to the the rainbow 
Schwinger-Dyson, and the need to perform a calculation 
in Minkowsky space is motivated. In Section ITTT1 the in- 
tegral in the mass gap equation is computed analytically 
in the case where the quark mass is assumed to be con- 
stant. The solution of the Minkowsky non-linear integral 
mass gap equation is addressed in Section llVl Finally In 
section the results and conclusion are presented. 



II. USING STANDARD APPROXIMATIONS 



and the mass gap equation in the instantaneous approx- 
imation is, 



M(p 2 



d|qlE" 

i 

M(q 2 ) 



4tt| P | 



log 



(IpI 



+ A 2 



UIpI 



+ A 



v/q 2 + M2(q2) 



(8) 



This is a one-dimensional non-linear integral equation. It 
has no singularities and it is solvable numerically Q. 

In the Euclidean approximation, it is assumed that the 
time component of the momentum can be replaced po — * 
ip4 or — po 2 — > p 2 . A convenient angular description of 
the variables is, 



In the literature the mass gap equation is usually 
solved either in equal time or in Euclidean space 
In the equal time approximation, the Lorentz invariance 
is lost. The space and time components of physical con- 
stants, say /„.* or / w * may differ [|J. Moreover it is not 
clear how to boost the hadrons outside the centre of mass 
frame. In the Euclidean approximation it is not clear if 
a simple Wick rotation p — > i Pi is exact because there 
may exist poles in the path ot the po axis. Moreover 
it is very hard to rotate back to Minkowsky with the 
inverse rotation i p^ — > po when only a numerical expres- 
sion of the dynamical mass is known. This is connected 
to the topological problem of discretising a curved sur- 
face. Although these two methods are approximate, they 
are used in the literature because they are fully consistent 
with dynamical symmetry breaking, providing the same 
approximation is also used in the bound state equations 
0| . At the same token these methods avoid the techni- 
cal problem of addressing poles and complex quantities, 
which are expected in Minkowsky space. By providing a 
solution of the mass gap equation, both Lorentz invari- 
ant and in Minkowsky space, both the equal time and 
the Euclidean time approximations may be better un- 
derstood. Here the approximate methods for solving the 
Schwinger-Dyson equation are reviewed. 

In the instantaneous or equal time approximation, it 
is assumed that the dependence in po and in go is ir- 
relevant in the kernel of the mass gap equation. The 
angular integral of the coulomb potential with neglected 
time component is, 



duj- 



2nq 2 ai 



i P' 
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27rq 2 ai 

z 2|pN 



log 



q 2 -2|p||q|w + A 2 
(|p|-|q|) 2 +Af- 



(|p| + |q|) 2 + A? 



(6) 



where the notation for quadri-vectors and for tri-vectors 
is p = (po,p). In what concerns the mass, a solution ex- 
ists where the mass is independent of go, and the integral 
in go is trivial, 



dq Q - 



<Zo z 



[q 2 + M 2 (q 2 ) 



^q 2 + Af 2 (q 2 ) 



(7) 
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(9) 



where the Euclidean momentum is, pe 2 = —p 2 = Pi 2 + 
P2 2 +P3 2 +P4 2 ■ In this case the boson exchange potential 
of eq. (0 is finite for all momenta. Covariance allows the 
choice of the external momentum parallel to the fourth 
axis, and three angular integrals can be performed, 



-AnoiqE 
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2ir aq — 



10 PE^ 
+ q 2 - A 2 + V(P 2 + Q 2 - A 2 ) 2 - Vg 2 



4p 2 g 2 



Ap 2 



^J-(p + q) 2 + \ 2 - ^-{p-q) 2 + \ 2 (10) 



which is a function of the real and positive variable — g 2 = 
qE 2 - Again the mass gap equation is reduced to a one- 
dimensional non-linear integral equation, 



M(p 2 



d(-q 2 )Y. 



2tt P 2 



p +q — X + 



P 



A 2 



Ap 2 q 2 



M{q 2 ) 
-q 2 + M 2 (q 2 ) 



(11) 



which has no singularities and is solvable numerically [j] . 

A third perspective, differing from the equal-time ap- 
proximation because it includes the retardation in the 
kernel and differing from the Euclidean approximation 
because it addresses time-like momenta, is provided by 
the coupled channel approach. Here the approximation 
consists in only considering the positive energy poles of 
the propagators. For simplicity I consider here a single 
Yukawa term, 
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Vq 2 + M 2 



P - q) 2 + A 2 



(12) 
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then the integral in go produces, 

dqo i i 



2n q - Hf + ie E — q — Hb + ie 



E- H f -H b + ie 



(13) 



The ttyj — . TTT3-, and the fermion mass M can be included 

Zn f ' Zri b ' 

in the boson creation vertex, 



rt 



M 



AH f H, 



(14) 



then the mass gap equation is equivalent to a coupled 
channel Hamiltonian equation, where the one fermion 
channel is coupled to the fermion plus boson channel by 
the boson creation and annihilation vertices, 



H f -ie-E T 

Tt H f + H b ~ie 



E 



cj)f 

<t>f,b 



0. (15) 



Reducing the equation by substitution of the fermion plus 
boson wave-function, the secular equation is, 







E-H f (p) 
1 



d 3 p 



(2tt) 3 

— rt(p, g)l 0(p)=O, (16) 



E-Hf(q)-H b (p-q) 



where it is clear that when E > M + Aj, the threshold 
for a boson production is open. This produces a cut 
in the function M(p), with an imaginary component for 
the mass. Actually more cuts appear when the energy is 
further increased, and the system couples to to 1 fermion 
and 2 bosons, 1 fermion and 3 bosons ... Therefore expect 
cusps and imaginary components are expected to appear 
when E = M + A;, E = M + 2\, ... 

When the mass gap equation is solved in the 
Minkowsky space, then all the features of the three differ- 
ent approaches of eqs. (JSJ, (fTTfl and (|Pp|l are expected to 
appear. Non-trivial solutions of the mass gap equation, 
and imaginary masses above thresholds are expected. 
Moreover It will be interesting to study the effect of the 
negative energy components on eigenvalues. 



III. ONE LOOP ANALYTICAL CALCULATION 

As a first approach to the Minkowsky integral, I as- 
sume a real constant mass M and compute analytically 
the integral of eq.JSJ). This can be regarded a one loop 
approximation to the mass gap equation, in the sense 
that the correct solution can be obtained iteratively, in 
an infinite loop calculation. The analytical result may 
also be eventually used to remove poles form the numer- 
ical iterative program. In this sense I choose to perform 
the multiple integrals in the same ordering that may be 
used in a numerical integration. A first attempt to start 



by an analytical integral of the three angular integrals 
(both trigonometric and hyperbolic) was abandoned be- 
cause the result was indeterminate. Therefore the two 
angular integrations of the three dimensional space are 
first performed, in order to simplify the kernel. Finally 
the double integration in the temporal go and in the spa- 
tial g are performed. 

The main aim of this section is the covariant study 
of the positive p 2 case which was not accessible to the 
approximate approaches. Covariance allows one to con- 
sider the p = (p°, 0) case. Analyticity can be used later, 
to continue this function to a p 2 negative case. The two 
dimensional angular integrals are trivial. The integral in 
go is also directly computed, using the residue theorem, 



J(P0,0) = J d 3 qd qo J2 



t (qo ~Po) 2 - q 2 - Ai 2 + ie 



tq* 



qo 2 - q 2 - M 2 + ie 



-irAirq 



2^q 2 + A " - * e Vq 2 + M 2 — ie 



^p + \Jq 2 + A 2 - ie + yjq 2 + M 2 - ie 
1 



=(}7) 



-po + y/q 2 + Xf-ie + y/q 2 + M 2 - ie 

The integral in q — |g| can also be performed analyti- 
cally. It is interesting to remark that the threshold for 
an imaginary contribution appears when pq > M + Aj. 
Indeed the pole appears at the root p of, 



Po 



= Vq 2 + m 2 



A 2 ^|q|-p, (18) 



P = 



'(p 2 -M 2 -X 2 ) -4M 2 A? 



4pt 



where p coincides with the mass shell momentum above 
threshold. A simple algebraic simplification transforms 
the integral into, 



d|q| 



M 



f2 



-7T 47rq 2 



%>o q - P I \JM 2 



Pi 



At 



-M 2 +p 2 + A 2 



(19) 

where the pole p defined in eq. (|19fl is real when the 
threshold opens, at p > M + A,, p is also real when 
Po < | — Ai|, below the pseudo-thershold, 



P 



b 2 -(M + A i ) 2 ][p 2 -(Af-A I )- 



ic 



(20) 



however in this case the residue vanishes, therefore I is 
only expected to possess an imaginary component above 
the threshold for the fermion-boson production. 

After cumbersome calculations, where the ultraviolet 
divergent terms cancel because they are proportional to 
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Qfj, the integral I(po, 0) can finally be reduced to the exact form, 
I 



Hp 2 ) = E 



2p 2 



9 9 ON M 2 

-M 2 +p 2 + A?)log^- 



-M 2 +p 2 + X 2 ) 



+(M 2 +p 2 -\f) 



.vVTa 2 
p 



log 



log 



log 



-log 



i 1 + 



vV + M 2 



(21) 



where p is defined in eq. i)19[l . The result of eq. I|21l) is 
not only correct for time-like momentum, it also applies 
to negative p 2 ,, where the integrals of eq. (|17|l remain 
correct. The integral I(p 2 ) is depicted in Fig. ^ 

There are three particular cases of eq. (I17fl that can be 
tested independently. The imaginary part can be directly 
computed with the residue theorem applied to ea. (|19ll . 
When pq > M + Xi the imaginary part is, 



Im[I(p ,0)}=J2- 



-Cti TT 



^{p 2 -M 2 -X 2 f-AM 2 X 2 



Pi 



Another important particular case is the matching 
point, between the time- like momentum and the space- 
like momentum, of p 2 = p 2 , = 0. The integral of eq. i|17[) 
is then easy to compute, 



7(0,0) = E~ a ^ 



2 M 2 log(Af 2 )-A 2 log(A 2 ) 
M 2 - A? 



(23) 



Finally in the space-like case p — (0, p) can be con- 
sidered. This implies that the poles are on the correct 
quadrants to enable a trivial Wick rotation. I can also 
use the angular integrals already performed in eq. 1)11) [I. 
The resulting integral can be computed analytically, 



J(0,P) 



2p 2 



(M 2 



A 2 ) log 



M 2 



(p 2 + M 2 + A 2 ) 2 - 4M 2 A 2 log 



(24) 



M 2 



p 2 + A 2 + V(M 2 + P 2 + A 2 ) 2 -4M 2 A 2 



M 2 + p 2 + A 2 - v /(Af 2 + p 2 + A 2 ) 2 -4M 2 A 2 
These three particular cases comply with eq. (|21|l . 



IV. APPROXIMATE ANALYTICAL SOLUTION 
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FIG. 1: Real and imaginary part of the integral I(p 2 ) that 
dynamically generates quark mass. Here the parameters M = 
0.3 GeV , X 1 = 0.1 GeV and A 2 = 0.2 GeV, on = 1 and a 2 = 
— 1 are used. The cusps occur precisely where the thresholds 
open, at p 2 = M + Ai and at p 2 = M + A 2 . 



Here the mass gap equation is solved in an analyti- 
cal one loop approximation. The standard method to 
solve the mass gap equation |J5J is the iterative method, 
where one starts by an initial educated guess for the mass 
Mi(p 2 ). I consider a constant Mi, which allows the an- 
alytical computation of the integral in eq. 0, with the 
techniques used in section II I II This produces the next 
mass in the iterative series, 



Mi 



M 2 {p 2 ) 



4-7T 



M 1 I(p 2 ) , 



(25) 



where the integral I{p 2 ) is defined in eq. 1)21(1. To find an 
exact solution to the mass gap equation, one would then 
need to continue this iterative process, computing again 
the integral with the function Mi{j> 2 ) — * Mz(p 2 ) and so 
on until the method converges. However these further 
iterative steps would probably need a numerical compu- 
tation, since the function M^p 2 ) is already a complicated 
one. Therefore, in this analytical approach, I choose to 
stop at the second step of the iteration. To minimise the 
error of this one loop computation, I demand that the 
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FIG. 2: The initial mass Mi, the one loop mass M2(p 2 ) and 
y/p 2 are illustrated. They all coincide at 0.3 GeV, and the 
parameters and ai = — = 36 are adjusted accordingly. 
The parameters are Ai=0.055 GeV and A 2 =0.105 GeV. This 
case also includes a second solution of the mass gap equation, 
M 2 = p 2 = (0.18 GeV) 2 . 



mass M 2 (p 2 ) coincides with the mass Mi at the mass 



shell momentum M 2 {p ) = p 



M 1 =M 2 (p 2 )=p 2 



(26) 



This is equivalent to assume that the mass dependence 
of the integral in eq. is dominated by the pole neigh- 
bourhood. The quark mass, dynamically generated in 
the mass gap equation, is expected to coincide with the 
constituent quark mass. The constituent quark mass is 
estimated in the quark model, where it is a crucial param- 
eter to produce the hadronic spectrum, and where it is of 
the order of 0.3 GeV. Therefore the coupling constants 
a, of the quark interaction J2J are adjusted to reproduce 
M2 = 0.3 GeV. An example of this is illustrated in Fig. 

I now discuss the different scenarios for the parame- 
ters Ai of the potential. I start with the case with three 
Yukawa terms. In the particular case of in eq. ©, this 
corresponds to the linear potential, which vanishes in 
the infrared limit. It occurs that in this case the gen- 
erated quark mass M 2 (p 2 ) is quite small close to the ori- 
gin p 2 ~ 0. It is then very difficult to arrive at a quark 
mass of 0.3 GeV. This may be related to the infrared 
cancellation. 



in 



d 3 k 



(2tt) 3 V k 2 + A? 



in 



in 



'V + X 2 



k 2 + A 2 



(27) 



which occurs when the Ai are equally spaced, A3 — A2 = 
A2 — Ai. Moreover the structure of the potential is quite 
complicated above threshold, with large cusps. Therefore 
this class of models is abandoned. 

Next the simpler case of two Yukawa terms is studied. 
This is related to the first derivative of eq. J3J. In this 
case it is easy to adjust the strength parameters ol\ = 



—a 2 to produce a quark mass of 0.3 GeV at the pole 
position of p 2 — 0.3 GeV 2 . This seems to agree with the 
infrared finite result, 



-in 



in 



d 3 k 



(2tt) 3 Vk 2 + A 2 k 2 + A 



Ai — A2 
in 



(28) 



Nevertheless, in this simple case of two Yukawa terms, 
there are two different scenarios. Depending on the steep- 
ness of the function M 2 (p 2 ), it intercepts the function p 2 
either at a single point or in two points. From a numeri- 
cal exploration one concludes that the two scenarios are 
separated by the line Ai(A2 — Ai) ~ 0.2 GeV 2 . This line 
is depicted in Fig. in the parameter space A2 Ai . There 
is a single solution when 



Ai(Aa-Ai) >0.2GeV 2 



(29) 



For instance this includes the case of a large A2, which 
acts as an ultraviolet Pauli-Villars cutoff, say Ai = 
0.3GeV, A2 = 3GeV. When Ai and A2 are quite large, 
the solution M(p 2 ) of the mass gap equation is very 
smooth, and the Instantaneous or the Euclidean ap- 
proaches constitute very good approximations to the ac- 
tual Minkowsky solution. It is also clear that in the case 
of a Pauli-Villars regularisation A ~ I is not acceptable. 
A would need to be computed, nevertheless its compu- 
tation could be performed with the techniques presented 
in this paper. 

One also concludes that there is a double solution when 



Ai(A 2 - Ai) < 0.2GeV 2 



(30) 



This includes the case of small and similar Ai. This case 
is related to the first derivative of eq. 10}. In the limit of 
vanishing parameters Ai, eq. 0} corresponds to a linear 
potential with a negative and infinite constant shift, 



err 



(31) 



3 

2.5 

2 

1 . 5 

1 

0.5 




ONE SOLUTION 



0.5 



1 . 5 



2 



2.5 



Ai 

~3~ 



FIG. 3: Phases with one and two solutions of the mass gap 
equation JSJ. 
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Comparing with the potential defined in eq. J3J the 
string tension is er = «i(A| — Af)/2. Assuming a string 
constant of the order of 0.14 GeV 2 estimated from the 
quark spectrum, one arrives at the following parameters, 



a(\ 2 - Ax) = 1.8 GeV , 
= 0.08 GeV , 



(32) 



and this implies that the potential also includes a nega- 
tive constant shift of the order of -1.7 GeV. An example 
of the generated masses is depicted in Fig. |5J 

To proceed iterating the equation (J3J), say at two loop 
order or more, or to compute the quark condensate {4>ip), 
the analytical integral becomes quite complicated. A 
numerical study is probably necessary and this is not 
addressed in this paper. Nevertheless the qualitative 
changes to the one loop computation can be anticipated. 
In what concerns the positive p 2 , it is expected that new 
channels will open whenever p 2 > (M + n\\ + mA2) 2 , 
where m and n are positive integers. This will affect 
the imaginary part of the mass. In what concerns the 
behaviour of the mass M(p 2 ) for very large positive or 
negative p 2 , an inspection of eq. J5J) shows that 



M(j?) 



An ai\2 — a^Ai 



(P 2 ) 2 
i d A q 



M(q 2 ) 



_ x (2tt)4 q*-M 2 (q 2 ) + ie 



(33) 



and therefore the mass should vanish proportionally to 
l/(p 2 ) 2 . In the one loop approximation the mass only 
vanishes like l/p 2 , see Fig^ Therefore the large momen- 
tum behaviour of the generated quark mass is expected 
to improve in the next iterations. 



component. In the case of a linear potential the mass 
gap equation is expected to have at least two solutions, 
and this agrees with references Q. 

The analytical continuation of a the numerical Eu- 
clidean space solution (with negative p 2 ) into the full 
Minkowsky has been studied in the literature Here 
I verify that this analytical continuation is not uniquely 
defined, even when a dense set of points is known in the 
the space-like p 2 < sector. In particular the equation 
(jSJ) includes at least three vanishing imaginary numbers 
— ze, summed respectively to the masses Ai, A2 and M. 
For external space-like momenta p 2 < 0, these ie are irrel- 
evant. However for time-like momenta p 2 > there are 
possible branch cuts both above threshold p 2 > (M+A;) 2 
and below the pseudo-threshold < p 2 < (M — A^) 2 , and 
the integral depends on the sign of each of the three van- 
ishing e. Only one integral, with all the three e > 
is causally correct. For instance in Fig. ^ the correct 
imaginary part of the mass is negative, but a continua- 
tion with a positive imaginary would also be analytically 
possible. As another example, a naive extension of the 
Euclidean eq. I|25|) to p 2 > fails to coincides with the 
full Minkowsky solution eq. 1)21(1 for some time-like mo- 
menta. 

The next step of this study will consist in applying nu- 
merical methods to proceed with the study of the mass 
gap equation. A numerical integration seems to be nec- 
essary to compute the quark condensate (i/jtp) ■ Moreover 
the exact numerical solution of the mass gap equation (JSJ) 
may also be attempted. With the analytical method of 
this paper, it may be possible to subtract the poles from 
the integrand and to compute their contribution analyt- 
ically. Another interesting method may be the fully nu- 
merical Monte-Carlo method [8j. These numerical meth- 
ods will be applied elsewhere. 



RESULTS AND CONCLUSION 



I address the technically challenging problem of solving 
the non-linear integral Schwinger-Dyson equation in full 
Minkowsky space. An analytical approach is followed, 
and approximate but analytical expressions for the quark 
mass are obtained. 

I find that the quark mass exhibits a branch cut above 
the threshold for boson creation, including an imaginary 
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